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Abstract: Adopting the horizon-crossing approximation, we derive the spectral index of 
fNL hi general N-flation model. Axion N-flation model is taken as a typical model for 
generating a large /jyx which characterizes the size of local form bispectrum. We find that 
its tilt rif NL is negligibly small when all inflatons have the same potential, but a negative 
detectable nj NL can be achieved in the axion N-flation with different decay constants for 
different inflatons. The measurement of nj NL can be used to support or falsify the axion 
N-flation in the near future. 
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1. Introduction 

In the single-component slow-roll inflation the different Fourier modes of the curvature 
perturbation are roughly uncorrelated and their distribution is almost Gaussian [|lj. Using 
8N formalism ^j, the curvature perturbation can be expanded to the non-linear orders 

( = 5N = Y^ NiSfa + \Y1 • v <;^ , >'; + lYl NijkWMjWk + ... , (1.1) 

2 . . D . . 

i 1,3 t,j,k 

where the subscript ,i denotes the derivative with respect to fa. The above formula indi- 
cates that the curvature perturbation is determined by the evolution of the homogeneous 
universe corresponding both to the classical inflation trajectory and to nearby trajectories. 
In principle, an important difference for the multi-field model from the single-component 
case is that the classical trajectory is not uniquely specified by the potential, but rather 
has to be given as a separate piece of information. It provides an opportunity to gen- 
erate a large non-Gaussianity in multi-field slow-roll inflation even at the inflationary 
epoch [||, |], ||, not at the end of inflation [||, ^, ||]. See some other relevant references 

in @ 0, 0, 0, m, @, m, u, m. 

On the other hand, not only the size of bispectrum f^L is measurable, but also its 
scale dependence characterized by n f NL is possibly measured by the accurate cosmological 
observations, such as Planck and CMBPol. In []T9f the authors concluded that Planck and 
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CMBPol are able to provide a \-a uncertainty on the spectral index of fNL for local form 
bispectrum as follows 



50 1 

A ri fNL ^0.1- == for Planck, (1.2) 

JNL yj J sky 

and 

50 1 

An fNL ~ 0.05- 7 = for CMBPol. (1.3) 

JNL \J fsky 

Recently the scale dependence of Jml for inflation and curvaton model is discussed in 



|20| , 21, 22, 24] and a large parameter space for generating a detectable nj NL in curvaton 



model is illustrated in [22, 23]. The value of nt NL will be an important discriminator to 
distinguish different models. 

In this paper we focus on a class of multi-field inflation, so-called N-flation, in which 
there are no cross couplings among the inflaton fields or the cross couplings are negligibly 
small. For an example, we consider the axion N-flation with a detectably large Jnl- But 
we find that nj NL is negligibly small when all inflatons have the same potential. However, 
a negative nj NL with detectably large absolute value can be obtained in the model with 
different decay constants for different inflatons. 

Our paper is organized as follows. In Sec. 2 we derive the spectral index of /nl for 
general N-flation. We explore the observables, including rif NL , in the axion N-flation in 
Sec. 3. Some discussions are contained in Sec. 4. 

2. Scale dependence of f NL in general N-flation 

In this section we will derive the spectral index of Jml in general N-flation with Nf un- 
coupled inflaton fields. The potential of inflatons is given by 

N f 

V = Y,Vi{<t>i). (2.1) 

i=i 

The dynamics of inflation is governed by 



N 



& = \i^(\$ + V i {<t> i )\ (2.2) 
+ m& = —V-, (2.3) 



where V- = dVi/dfa. In this paper we work on the unit of M p = 1. In the slow-roll limit, 
the above equations are simplified to be 

H 2 * j, (2.4) 
3ff& ~ -V(. (2.5) 
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Similarly, we introduce some 'slow-roll' parameters, such as 

i fvr 2 



From Eq.(|2.4|), we find 



and then 



or equivalently, 



^2 UH • (2 ' 6) 



12 

i ' i 



H 2 2 V 2 ' v ' 

en = £(|) 2 e, (2. 



Inflation happens when en < 1. The multi- field slow-roll inflation can be achieved even 
when some inflaton fields have steep potentials because the friction term for each inflaton 
is contributed from all of fields. 

Following [Pq[ , the number of e- folds before the end of inflation is given by 

■c-L fPi V- 

N = J2 rT/dfa. (2.10) 

~T J(h end v i 



t=l " ^ 



Therefore we obtain 



and 



^i=(l--^-)%- (2-12) 



Using the SN formalism |2|, the power spectrum of scalar perturbation is given by 

* = (£f)"E*3- < 2 - 13 > 

i 

A slow-roll inflation model predicts a near scale-invariant power spectrum of scalar per- 
turbation. The deviation from exact scale invariance is described by the spectral index n s 
which is defined by 



dlnP C . . .,>:, , '-,-V,,V. ; 



ns = 1 + : tE^ = 1 - 2e H - " (2.14) 



where 
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See [25, 26] in detail. In this paper, we assume, without loss of generality, cfii < 0, so that 
V- > for V i. On the other hand, the gravitational wave perturbation is also generated 
during inflation and the amplitude of its power spectrum is determined by the inflation 
scale 

ft = ^, (2.16) 
and then the tensor-scalar ratio becomes 

r T = P T /P c = 16 /V- • (2.17) 

/ i 6i 

The amplitude of gravitational perturbation is also near scale-invariant and its tilt is defined 
by 

d\nP T 

n T = -j^j- = ~2e H . (2.18) 

In single-field inflation, r-r = 16e# which leads to a consistency relation ut = — ry/8. For 
iV-field inflation, it is modified to be 

n T < -r T /8. (2.19) 

Here we used the inequality 

e H > 1 /J2 ~ ' ( 2 - 2 °) 

/ • ei 

where the equality is satisfied when €{ = ejfV/Vi for V i. 

From Eq.tfLll), the non-Gaussianity parameter /jvl can be written by 

fm= (2 ' 21) 

Since iVy oc 6™ for N-flation, the above formula is simplified to be 

fNL = J2f l NL, (2-22) 



where 



_B NjN u _5r| 1 / o 
^ L ~ 6 (J2 k N\) 2 ~ 6 128 ^ V '• 



Following [21], the spectral index of f^L for N-flation is 

Inl 



a f x , -y—J2^NL( 2n muiti,i + n f)i ), (2.24) 
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where 



9 >r- ( N k Nj N k \ 

n mu lti,i = 2 Z^/ ekl \ v AT 2 ~ il W ' {2,. lb) 

k,l \^i iV ,i >V 
AT i 



k 

and 



f2> oVl V" VL V'V" V'" 

F, (2) = -2(^) 2 -± + + 2-±-^5 ki - ^5 ki . (2.27) 

kit K V V VV VV V 

Considering iVj = Vi/V^, we obtain 



/>,„„/„., (l- "^p- ) , (2-28) 



Using Eqs.(U|), (fT23|) , ( gjgg ), (|2~28|) and ( gjgj ), we can easily calculate the spectral index 



of /atl in general N-flation. The first term in ( 2.28| ) contribute 2(1 — n s — 2en) to nf NL . If 
e H < 1, 2(1 - n, - 2e H ) ^ 0.08 for n s = 0.96. 

In addition, we are also interested in the trispectrum in this case. See the Appendix. [A| 
in detail. 

3. Axion N-flation 

In this section we focus on axion N-flation model || in which there are Nj inflaton fields 
and the potential of fa is given by 

Vi = Af(l-cosQi), (3.1) 

where Oj = 2ir4>i/ fi and /, is the i-th axion decay constant. The 'slow-roll' parameters e^'s 
are given by 

_ 2tt 2 1 + cos otj 

ff 1 - COS Oii 



and 



VV" 1 

N ,a = l -^7t = ^ • ( 3 - 3 ) 



Vf 2 1 + cos oti 
Now the number of e-folds before the end of inflation becomes 



N ~y(l£\ i n 1 . (3.4) 

^>\2ir 1 + cosaj v ; 

4 = 1 
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Quoting the results in ||, the observables of axion N-flation are 



8vr 2 v-^ Af 1 /v-^ 1 , . 

r T = 16/J2-, (3-7) 

. 5 rf, ^ 1 1 

6 128 f-i 1 + cos Oj 

Using the results in the previous section, we can calculate the spectral index of f^L for the 
axion N-flation, namely 

n mu iu,i = (1 - n s - 2e H ) - ^pj$, (3-9) 

n u = <yfti- 2 y^ (3-10) 



and then 



, 16vr 2 ^ Af 1 1 /r-y 1 1 

n /w£ = 2(1 - n s - &*) -3^^ ^^l + CO sa 4 / X> 7, TT^osaj 



.1 



+ 2 y( 2 (^) 2 -^ — l - — /yl — 1 — . (3.11) 

^\ K V> VJ l + cosai/ ^ ej 1 + cosaj- V 7 



The first and second lines in Eq.( |3.11 ) are contributed by n mu i ti ^ and n^, respectively. 



Since the second term in n mu \ti,i is large compared to n^j for ~ 7r which corresponds to 
a large /jvz,, the spectral index of /nl is mainly determined by n mu iu t i. 

3.1 Axion N-flation with the same potential for different inflatons 

In this subsection we discuss the model with Aj = A and fi = f for V i in detail. Now the 
tensor-scalar ratio and /nl become 

rT = — /iZTT^r (3 - 12) 

107T 2 ^ 1-COSOi / |^ 1-COSOj I 

" "37T (1 + cosaj) 2 / 1 + cosa, 1 ' ( - K l0) 
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and the spectral indices of scalar power spectrum and /jvx are simplified to be 

8^ 2 A 4 , N 

n s = l-2e H - W tp, (3-14) 



I 2 

on - 

i 

5rx 1 



/^(l-n s -2e // ) 2 ^l(l-cosa i ) / ^ -— 

8ir z ' Si / ' e-j 1 + cosa,- 

i I 3 

(l-n s -2e H ), (3.15) 



96 f NL 



where 



A 4 / 

and 

2"7T 2 v— v 1 — COS 2a; / I v-^ I 

^ = JtH 2 l / -<*»<*,)] • ( 3 - 17 ) 

We see that fimuiu i 

= and rif NL is expected to be small in this simple setup. We will show 
the smallness of nj NL and extend the discussions in || in the following two subsubsections. 



3.1.1 AL ~ A, 



tot 



Assuming the initial conditions for a% satisfy a uniform distribution, the number of e-folds 
before the end of inflation is roughly given by 

Ntot * 7^f N f- (3-18) 

A large enough total number of e-folds, for example not less than 60, can be achieved as 
long as the number of inflatons is large enoug h, namely Nf > 1.7 x 10 3 // 2 - 

Similar to ]|], we consider the case with A* ~ Ntot in this subsubsection, where AT* 
denotes the number of e-folds corresponding to the CMB scale. So we have 

N ,f^ (, 19 > 

and A 4 /(3-£f 2 ) ~ 1/Afy, (h — 7I " 2 /(A^// 2 ). The spectral index of the scale power spectrum 
is related to A* by 

n s * 1 - (3.20) 

which is the same as that in ||. For AT* — 60, n s ~ 0.942. The numerical results for 
A* ~ Ntot — 60 and / = 1/2 are illustrated in Fig. [|. The spectral index of scalar power 
spectrum is nicely compatible with the estimation, and the tensor-scalar ratio is roughly 
0(1O~ 5 ~ 10~ 3 ). As long as Nt is large enough, from Eqs. ( p. 12 ) and (3.13), both vt and 
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io- 13 



10 20 30 40 50 60 

/nl 

Figure 1: The observables of axion N- nation with the same potential for different inflatons. Here 
/ = 1/2 and the initial conditions of on is assumed to be the uniform distribution. 

/nl are proportional to l// 2 . From the numerical results, we find that f^L is bounded 
from above, namely 

57T 2 

fNL < ^ ' ( 3 - 21 ) 

For / < 1, the value of f^L can be of order O(10) which is large enough to be detected in 
the near future, but its tilt rif NL (— O(10 -7 ) ~ — 0(1O -14 )) is too small to be detected in 
this case. 
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3.1.2 A few inflatons stays around the hilltop 

In ||, the authors proposed an alternative way to achieve a large /nl- Roughly speaking, 
the summation of the function which is proportional to 1/ej is dominated by those fields 
with the smallest ej, or equivalently staying around the hilltop. Suppose some number N 
of fields have roughly comparable ej, of order e which is related to the angle a by 

IT 2 , , 7T 2 



e~^(l + cosa)~^(7r-a) 2 , (3.22) 
with (tt — a) -C 1. In this case, P^, rx, and /nl becomes 



H 2 N H 2 Nf 2 
P< * T-o^ * 7-17— 4^ (3-23) 



2 I 4vr 4 (7r-a) 2 ' 



Vol 8vr 2 . n2 
r T ~ ~ — ^-(vr - a) 2 , (3.24) 
iV Nf 2 

/iVL - (3.25) 



Considering N > 1, 



57T 2 

/jvl < ^ , (3.26) 

which is the same as that in Sec. |3.1.1| . Now the spectral index of /jyx is given by 

1 + cos a . , ,„ .. 

r 2 (1 ~n s - 2e H ), (3.27) 

which is much smaller than (1 — n s — leu). So the spectral index of /nl is undetectable 
in this case as well. 

Because the potential become very flat around the hilltop, one may worry that the 
classical motions of those inflatons staying around the hilltop are not reliable. Requiring 
that the classical motion of inflaton (f> = fa/2ir per Hubble time be not less than its 
quantum fluctuation H/2-k yields 

(7r "" ) ~iS = 2rl 7' (3 - 28) 

where £ = 8vr 2 A 4 /(3i7 2 / 2 ) = (l-n s -2e H ) < 1. 1 On the other hand, from ( pT23| ), (vr-a) 
can be estimated by 



where An = \/P£. Therefore 



(vr-«) = ^, (3.29) 



V^f > (3.30) 



1 In jH[, the authors consider \4>huito P ~ <t>\ }t H/2ir, namely (it — a) > H/f, which is much looser than 
what we require, where <j>hUltop = //2. 
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Taking Eq.( 3,25| ) into account, /nl is bounded from above by 

5 £ LfTr 



Inl < 



12 A n 

2 _ o An-, in— 9 



(3.31) 



WMAP normalization Q implies P c = = 2.461 x 10~ 9 , and then f NL < 8 x 10^/VN. 
Since £ < 0(1CP 2 ), /nl ^ 80 /y/~N. The value of /nl cannot be very large in this case. 

3.2 Axion N-flation with different potential for different inflatons 

In this subsection, we consider a class of axion N-flation in which Aj = A for V i, but the 
decay constants can be different for different inflatons. For simplicity, we only focus on the 
cases with only two different decay constants. 

3.2.1 One inflaton with different decay constant stays around the hilltop 

In this subsubsection we consider a model with potentials 



V{<j>i) = A 4 (l-cosai) 
V{4>) = A 4 (l - cos a), 



(3.32) 
(3.33) 



where aj = 2ir 4>i/f for i = 1,2,..., Nt and a = 2ircj)/f. Here / is assumed to be different 
from /. We consider Nf ^> 1 and the initial values of aj satisfy a uniform distribution. 2 
The total number of e-folds before the end of inflation is given by 



Ni 



tul 



ln2 ~ f 2 , 



2vr 



4:ir 2 1 + cos a 



(3.34) 



We assume <p stays around the hilltop of its potential. Requiring that the classical motion 
of inflaton cb is reliable leads to 



v 71 " amt > ~ 47t2A 4 : 



(3.35) 



where aim denotes the initial value of a. 

In this model, the Hubble parameter is roughly related to A by H 2 ~ NjA 4 /3. Because 
the inflaton staying around the hilltop almost does not move, we also have 



en 



In 2 



N 



tot 



P , 2 
— o l n r 

2ir 2 it — a 



The spectral index of power spectrum and f^L become 



In 2 



N 



tot 2vr 

81n2 



^ln^ 



1 + 4- 



n 



fNL 



N 



tot 



, 7 In ■ 



Cr(N f )N f + (f/f) 2 J ^ j 
Cr(N f )N f + ^ d r (N f )N f + ^ 



(3.36) 



(3.37) 



c r (N f )N f + {f/fYj^ 



d r (N f )N f + (f/f)^ 



(it— a) 4 



(3.38) 



2 Thank C. T. Byrnes for the helpful discussion. 



where 



, » T , 1 v— v 1 — COS OLi , , „ rt 1 1 — COS OLi 

c r (N) = — > -, and dr(N) = — > - ^, (3.39 

v 1 A^^l + cosai' v ; TV f-' (1 + cosaj) 2 v y 

which depend on the total number of points and the random choice of the series V for 
CKj G [0,7r]. 

Similar to Sec. |3.1,1| , we consider that the total number of e-folds corresponds to the 
CMB scale, namely ~ Ntot — 60. Considering / < 1, the total number of e-fold is 
mainly contributed by 4>i for i = 1,2, ...,Nf. If / = /, n s — 1 — 51n2/iV* and ra^ Ari ~ 
which are consistent with the results in Sec. 



For / < /, the second term in the bracket of Eq. ( |3,37| ) is larger than 4 and then 
n s < 1 — 51n2/iV*. For / > /, we find ra s > 1 — 51n2/iV*. In general, it is difficult 
to get the analytical estimation for the model with f ^ f, but we can expect that the 
results should be quite different the model with the same decay constants. The numerical 
results are showed in Fig. 0. Here we keep / = 1/2 and (it — ai n i )/?r = 2 x 1CT 4 fixed and 
explore the results for two choices of /. From Fig. 0, we see that a detectable negative 



rif NL is possibly generated for / = 1/4. We also check the requirement in ( [3.35D which is 
satisfied as long as (it — a) /it > 2 x 10~ 5 for / = 1/4. Unfortunately, the spectral index 
n s (< 0.91) in the case with detectable rif NL is small compared to WMAP 7yr data |^7j 
(n s = 0.968 db 0.012). See the green dots in Fig. ^. For / = 1, the spectral index can fit 
the WMAP 7yr data nicely, but the scale dependence of f^L becomes negligibly small. 
See the red dots in Fig. ^. This setup indicates that a large rif NL can be achieved in axion 
N-flation. 

3.2.2 A general axion N-flation with two kinds of inflatons 

In this subsubsection we consider a more general axion N-flation model in which there are 
two kinds of inflatons who have different decay constants, 

Vi(<k,i) = A 4 (l-cosai,i), (3.40) 
V 2 (4>2, 3 ) = A 4 (l-cosa 2ii ), (3.41) 

where = 2ir(j>i i i/ f\ for i = 1,2, ...,Ni, and «2,j = 27r^2j//2 for j = 1,2, N 2 . Here f\ 
is assumed to be different from f 2 - We also consider N±, N 2 S> 1 and the initial values of 
an and a 2 j satisfy a uniform distribution respectively. In this model the total number of 
e-folds before the end of inflation is given by 



In 2 

2^' 



Hot^— (ffN 1 + flN 2 ). (3.42) 



In this model, the Hubble parameter is related to A by H 2 = (Ni + N 2 )A 4 /3 and the 
slow-roll parameter is given by 
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Log 10 r T 





Figure 2: The observables of axion N- nation with different decay constants. Here we adopt / = 1/2 
and (tt — ai n i)/^ — 2 x 10 -4 . The blue and green dots correspond to / = 1/4 and the red dots 
corresponds to / = 1. The green dots illustrate the case with detectable scale dependence of Jnl- 



Similarly, we can also derive the spectral index of power spectrum and /nl , namely 



27T* 



i>0 ~T~ nO 



57T 



c^N^ + c r2 (N 2 )N 2 



(Nt + N 2 y V fi /I J N 1+ N 2 c ri (NjffNi + c r2 (N 2 )f*N 2 ' 



16vr 2 



-V'l 



(N 1 )N 1 + c r2 (N 2 )N 2 



d ri (Ni)N 1 + dr 2 (N 2 )N2 



c ri (NjftNi + Cr2 (N 2 )f$N 2 d ri (iVO/fiVi + d r2 (N 2 )f$N 2 



(3.44) 
.(3.45) 
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Again, if f\ = /2, the scale dependence of /nl is negligibly small and our results recover 
those in Sec. |3TT| . 

For the case with f\ ^ /2, we need to adopt numerical method to explore it. Con- 
sidering iV* ~ N to t — 60, our numerical results are illustrated in Fig. 0. Here we adopt 



Log,o r T 
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X i • •• 
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' detectable region [CMBPol] 
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/nl 
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Figure 3: The observables of a general axion N-flation with different decay constants. Here we 
adopt /i = 1/2 and /a = 1/4. The green dots illustrate the case with detectable scale dependence 
of /nl- 
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fx = 1/2, JVi = 4800, f 2 = 1/4 and N 2 = 8000. We see that a detectable scale dependence 
of Jnl can be obtained. 



4. Discussions 



In this paper we work out the spectral index of /at l m general N-flation model by adopting 
the horizon-crossing approximation. For an instance, we focus on the axion N-flation 
proposed in || and we find that Jnl can be large enough to be detected in the near future. 

• All inflatons have the same potential. In the case with the total number of e-folds 
corresponding to the CMB scale, f^L is roughly of order O(10) for a uniform distribution 
of initial angles if the decay constant / is around the Planck scale. Our numerical results 
indicate that Jnl is bounded from above by 57r 2 /3/ 2 , and tnl — (2tt//) 2 ^ 3 /^x- Another 
way to obtain a large /jvl is that a few inflatons stay around the hilltop and the summations 
in the formula of /nl is mainly contributed by these inflatons. Typically Jnl is less than 
one hundred in the second case. In both cases the tilt of Jnl is negligibly small. The 
smallness of n j NL in these two case comes from the accident cancellation in n mu iti,i because 
all inflatons have the same potential. 

• Different inflatons have different decay constants. In this kind of model, the accident 
cancellation in n mu iti,i does not happen and a negative detectable nj NL can be obtained. It 
implies that a scale independent Jnl is not generic prediction of N-flation. A more general 
axion N-flation is left to be studied in the future. 

In addition, we also study the trispectrum for N-flation in the Appendix. |A|. We find 
that there is a universal relation between g^L and tnl, namely g^L — %5t~nl/27, for the 
cases with large Jnl- But the relation between f^L an d rjyi depends on detail of the 
model. 

If there is only one inflaton field, the axion N-flation becomes natural inflation. The 
slow-roll parameters are roughly given by M 2 // 2 and then the slow-roll conditions are 
satisfied if / > M p . However such a large decay constant cannot be achieved in some 
fundamental theories, such as string theory |28], |29|]. In [|30| the N-flation is supposed to 
realize a slow-roll inflation even when / < M p by calling multi inflaton fields. On the 
other hand, in [31] we proposed a conjecture that the effective gravity scale should be 
Ac = M p / y/Nf for the case with multi scalar fields, not M p . If the decay constant of 
axion field is required to be smaller than Aq, namely Nff 2 < 1, the total number of e-folds 
is less than one and the axion N-flation fails. 
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A. The trispectrum in N-flation 

From Eq. ( |1 . 1|) , the size of trispectrum is characterized by two parameters tjvl and g^L 
which are defined by 

TNL = J2 N ,ij N ,i* N ,j N ,k / (E^') 3 ' ( A>1 ) 
i,j,k I I 

9nl = § E KijkNiNjNk / C£ N J) 3 • ( A - 2 ) 

i,j,k I I 

For N-flation, the above two equations are simplified to be 

i I I 

= § E A -- A ' / (E A ') 3 - ( A - 4 ) 



where 

,y,y," 2 y» v^" 

V' 3 V V' 2 



= 2-i-±- (A.5) 



i 



We remind the readers that the value of tjvl is bounded from below by (^/nl) 2 - For the 
axion N-flation, we have 



2tt sin a,- 




(A.6) 



Therefore 

r - V -I' 1 ~ r "" n/ /|V- - 1 I (A 7 s ! 

'■ SL ^ l^il+muvi'i / I ^47r 2 l + cosaj I ' 

_ 25 y jj I 1 - cosa 2 ./ f V f f 1 -cosaA 
' x L " 54 ^ 4vr 2 (1 + cos a*) 3 / [ 4^ 4vr 2 1 + cos aj I ' { ' 

For the case with fi = f for V i, both t^vl and 5atl are proportional to l// 4 . Similar 



to Sec. |3.1.1| , we consider N* ~ AT tot ~ 60 and the numerical results for / = 1/2 are showed 
in Fig. |j. Fitting the numerical results, we find that t^l and g^L are roughly related to 
f NL and / by 

r 7VL — ( ~T" ] fNL> 9NL ^^T NL . (A.9) 

Roughly speaking, f^L ~ I// 2 and then g^L ~ §|t7vl ~ I// 4 which is consistent with 



our estimation. On the other hand, combining ( A.9 ) and t^vl > (§/tvl) 2 j we get Jnl ^5 



4tt 2 /3/ 2 which is roughly the same as our results in Sec. [3.1.1| . 
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Figure 4: The non-Gaussianity parameters of axion N-flation for the uniform distribution of on 
when / = 1/2. Here we consider the case with iV* ~ N to t — 60. 



The non-Gaussianity parameters of trispectrum for the case in which a few innatons 
stay around the hilltop were worked out in ||: t^l = (§/tvl) 2 and g^L = ^tnl, which 
are different from the previous case unless /nl — 47r 2 /3/ 2 . 

The non-Gaussianity parameters for the model in Sec. 3.2.1 are showed in Fig. || 
Fitting the green dots corresponding to a detectable scale dependence of /atl, we find 

25 

t nl ~ 7.78/^F, g NL ~ — tjvl. (A.10) 

Fitting the blue and red dots respectively, we reach the same results 

t nl ~ 8.28/^| 9 , g NL ~ frjvi. (A.ll) 

Here the relations between r^r^ and /jvl depend on the choice of the parameters, such as 
/, / and a, in the model. 

Similarly we also investigate the non-Gaussianity parameters for the model in Sec. 3.2.2. 
See the numerical results in Fig. |6[ Fitting the blue and green dots respectively, we find 

25 

t nl ^ 7-Wnl , 9nl ^ ^r NL . (A. 12) 

To summarize, there is a universal relation between g^L and t^l-, namely g^L — 
YfT^L- The value of g^L is the same order of t^l- But the relations between t^l and 
Jnl are different for different cases. 
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Figure 5: The non-Gaussianity parameters in the axion N-flation where Nf infiatons have the 
same decay constant /, and the decay constant of another infiaton is /. Here we consider the case 
with iV* ~ Ntot — 60 and adopt / = 1/2 and (ir — otini)/^ = 2 x 10 ~ 4 . The blue and green dots 
correspond to / = 1/4 and the red dots corresponds to / = 1. The green dots illustrate the case 
with detectable scale dependence of Jnl- 
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